Insensitivity of superconductivity to disorder in the cuprates 
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Using a dynamical cluster quantum Monte Carlo approximation, we investigate the effect of local 
disorder on the stability of d-wave superconductivity including the effect of electronic correlations 
in both particle-particle and particle-hole channels. With increasing impurity potential, we find 
an initial rise of the critical temperature due to an enhancement of antiferromagnetic spin correla- 
tions, followed by a decrease of T c due to scattering from impurity-induced moments and ordinary 
pairbreaking. We discuss the weak initial dependence of T c on impurity concentration found in 
comparison to experiments on cuprates. 



PACS numbers: 74.62.Dh,74.25.Ha, 74.72.-h, 74.81.-g 

Introduction. Disorder is an essential feature of the su- 
perconducting cuprates. Crystal growth procedures lead 
generically to defects such as grain boundaries, atomic 
site switching, and vacancies. Additional disorder, often 
in the form of oxygen or other charged defects, is almost 
always introduced away from the Cu02 plane upon dop- 
ing the parent compound from the Mott insulating state. 
This last type of disorder may be responsible for local 
nanoscale electronic inhomogeneity in the superconduct- 
ing state of the cuprate Bi-2212 indicated by scanning 
tunnelling spectroscopy (STS) experiments [l|, [2j, y, |j|. 
These experiments show modulations of the local gap 
near the impurity sites on the order of the correlation 
length [a, |aj. A recent experiment imaging high-energy 
resonances thought to be the dopant atoms themselves 
shows a strong positive correlation of the magnitude of 
the local spectral gap with the locations of the dopants, 
leading to suggestions that the origin of the observed gap 
modulations are caused by atomic scale variations in the 
pairing interactional, [6|. 

When impurities like Zn and Ni are substituted for Cu 
in the Cu02 plane, or planar defects created byelectron 
irradiation, superconductivity is suppressed [71. la. |{j Ilfj|. 
Because the screened Coulomb potential due to these 
defects is very short range llj , such impurities are fre- 
quently modeled by pointlike ((^-function) scatterers. 
The expected form for the suppression of superconduc- 
tivity in the BCS theory of d-wave superconductors is 
then identical to the expression given by Abrikosov and 
Gor'kov[12| for magnetic impurities in s-wave supercon- 
ductors (see, e.g., Refll3Th However, experimentally a 
significantly slower initial slope of the T c suppression is 
observed. For example, Tolpygo et alflfjl reported a sup- 
pression 2-3 smaller than the AG curve [10(| in resistivity 
measurements of YBCO films. Other unusual deviations 
from AG behavior have been observed at larger disorder 
levels; for example, an electron irradiation study[14[ on 
optimally doped YBCO reported a linear behavior in T c 
vs. resistivity over the entire T c range. 



Theoretically, several possible effects beyond 
Abrikosov-Gor'kov (AG) theory have been explored. 
A numerical mean field study of disordered d-wave 
superconductors [15j including the self-consistent sup- 
pression of the order parameter around each impurity site 
showed deviations from the AG result. Several authors 
attempted to account for the slowness of the T c suppres- 
sion by assunring that the scattering potential of planar 
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impurities was extended, or anisotropic [It 
Recently, Graser et al.[20( calculated both T c and the 
impurity resistivity p within a consistent model of 
extended potential scatterers, and concluded that the 
unusual T c vs. p behavior seen in cuprate experiments 
should be attributed to strong correlations or strong 
coupling corrections to BCS theory In general, the 
effect of correlations on the structure and scattering of 
quasiparticle states in a disordered d-wave superconduc- 
tors is still an open and very important question for 
cuprates and other unconventional superconductors. 

One interesting consequence of disorder in a correlated 
electron host is impurity-induced magnetism: nuclear 
magnetic resonance measurements indicate the forma- 
tion of magnetic moments upon chemical substitution of 
a nonmagnetic impurity for a Cu 2l|, |22| . This was cor- 



roborated by calculations of the magnetic spin suscepti- 
bility, which displays Curie- Weiss behavior upon impu- 
rity doping (see e.g. Ref. l23l ). Several aspects of the- 



ory and experiment in connection with disorder-induced 
magnetism in cuprates and ID spin systems have recently 
been reviewed in Ref. 24 . While most of the theoretical 



work on these questions has been confined to the normal 
state, the quasiparticles deep in the d-wave supercon- 
ducting state are also affected. Mean-field calculations 
utilizing the Gutzwiller approximation [25| suggest that 
the effects of disorder on the density of states are sup- 
pressed in the presence of strong correlations, specifically 
near the nodes and at low energies. Similar effects in the 
density of states are also recovered in calculations where 
correlations are treated in a simple Hartrce-Fock scheme 



by Andersen et al. [26| , who found however that although 
the effects of disorder on the density of states were indeed 
weakened, some unusual effects outside the framework of 
BCS theory were also present, e.g. the breakdown of 
universal transport in d-wave superconductors [27(. 

In this paper, we aim to understand some of the ef- 
fects of disorder on the suppression of the transition to 
d-wave superconductivity. First, a small concentration 
of weak impurities is shown to cause an increase in the 
effective antiferromagnetic exchange coupling, which en- 
hances superconductivity within the Hubbard model. At 
the same time, the disorder causes pairbreaking, which 
tends to suppress T c . As the impurity potential is in- 
creased, the pairbreaking overcomes the enhancement of 
J causing a decrease in T c , which continues until it sat- 
urates when the unitary limit is achieved. We suggest 
that these effects may partially account for the observed 
slow suppression of T c by disorder in the cuprates. 

Formalism. The Hamiltonian of our model is 

H = -t y^ cl a Cj C + U*^2 n i<? n i-<? + X! Vi n i<y (!) 

(ij)a i i 

where cT^cjo-) creates (destroys) an electron with spin 
a at site i, and rii a = c ia Ci a . Here (ij) denotes nearest 
neighbor sites i and j, U denotes the on-site Coulomb 
repulsion and and t is the nearest-neighbor hopping am- 
plitude. The impurity is modeled as a potential Vj = V 
on a single site. We shall give a brief review of the method 
(see Ref. [23| for further details). 

To study {T]), we employ the dynamical cluster approx- 
imation (DCA)[H, |2a, [30|]. The DC A is a dynamical 
mean-held theory which sclf-consistcntly calculates the 
self-energy on a cluster of size N c embedded in a host. 
Correlations on the cluster are treated explicitly. Interac- 
tions beyond the cluster scale are dealt with on a mean- 
held level within the self-consistent host. With increasing 
cluster size, the DCA systematically interpolates between 
the single-site dynamical mean held (DMFT) [3l| and the 
exact result. Cluster dynamical mean field calculations 
(including the DCA) of the Hubbard model are found to 
correctly obtain many of the features of the cuprates, 
including a Mott gap and strong AF correlations, d- 
wave superconductivity and pseudogap behavior [30]. To 
solve the cluster problem, we use a quantum Monte Carlo 
(QMC) algorithm [32], and employ the maximum entropy 
method [331] to calculate the real frequency dynamic spin 
susceptibility. The sign problem in QMC is small for the 
values of U considered, and is therefore not an issue for 
the calculations presented here. 

The result of the QMC calculation depends on the dis- 
order configuration, but T c is determined by the aver- 
age Green's function, which we compute in the following 
way. For a concentration x, contributions with m impuri- 
ties on the cluster are weighted by a combinatoric factor 
x m (l — x) Na ~ m . It is reasonable, for small concentra- 
tions (x < 1/N C ), to consider only those configurations 



with zero or one impurities. For the zero and one im- 
purity case, the combinatoric factors expand to I — xN c 
and xN c , respectively. We can then write the disorder 
average 



G% = xN c G c hij + (I - xN c )Gfa 



(2) 



where GL „■_• is the cluster real space Green's function for 
m impurities. The disorder-averaged Green's function is 
then used to continue the DCA algorithm. 

To determine the critical temperature T c , we extrapo- 
late the pair-field susceptibility Xd(^)|20]j an d n °te that 
the system enters the superconducting state when Xd(T) 
diverges. To interpret the results we present below, we 
will also need to calculate the induced magnetic moment 
rn. This is done using a method introduced by Krishna- 
murthy et al[34(. We note that the square magnetic mo- 
ment in the low-temperature limit is proportional to T 
times the magnetic susceptibility. To study the effect 
of the impurity, we subtract the pure susceptibility, and 
arrive at 



m induced <* MXl — Xo) 



(3) 



where Xi an d Xo are the susceptibilities of a cluster with a 
single impurity and a homogeneous cluster, respectively. 

Results. We carry out DCA/QMC calculations using 
the N c = 16, type A|35[ cluster for the Hamiltonian in 
Eq. [Q fix the doping at 10% and let U = At. Esti- 
mates for T c have been shown to be robust against clus- 
ter size effects 36[. Furthermore, we have investigated a 
possible finite size effect by observing the change zero fre- 
quency spin-spin correlation function (not shown) , which 
was found not to deviate appreciably from the clean clus- 
ter beyond the first nearest neighbour — indicating that 
the finite size effect does not play a significant role on the 
quantities we report. 

We first investigate the d-wave superconducting transi- 
tion temperature T c and the induced moment of the sys- 
tem as a function of impurity potential V for various 
values of the impurity concentration. 

Our first significant finding is the initial weak increase 
of T c in the region < V < t for 3% impurity concen- 
tration (Fig. [I]). This is completely unexpected from 
the point of view of AG theory, where any concentration 
of impurities of any strength will suppress T c initially. 
The increase in T c with respect to the homogeneous sys- 
tem is slightly less than 4%. After we increase V to 
a significantly larger value, for example 20£, the d-wave 
superconductivity still survives and the critical temper- 
ature saturates. This is consistent with the BCS theory 
of pair breaking by point like impurities of a d-wave su- 
perconductor (without correlations in the particle-hole 
channel), where increasing impurity potential past the 
bandwidth (~ it) drives the impurity into the unitarity 
limit where the scattering rate saturates [37[. Increasing 
the impurity concentration beyond 3% causes a dramatic 
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Figure 1: The critical temperature T c as a function of impu- 
rity potential for N c = 16 and U — 4£, at impurity concen- 
trations x = 3% and x = 6%. Error bars are calculated from 
the extrapolation of the pair-field susceptibility [231] . Inset: 
Blowup of the region of small impurity potential. 



monotonic drop in T c for all V > 0. For 6% impurity con- 
centration T r vanishes even before V = It. 
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Figure 2: The critical temperature T c as a function of impu- 
rity concentration for N c — 16, and U = 4i, at impurity po- 
tentials V — t, V = At and V = 20t. Error bars are calculated 
from the extrapolation of the pair-field susceptibility [231] . The 
AG result is a fit to the critical concentration for V = 20£. 

Fig. [5] shows the behavior of T c as a function of impu- 
rity concentration. For all concentrations considered, the 
initial slope is either positive or nearly zero, in marked 
contrast to the Abrikosov-Gor'kov curve, which has a 
negative initial slope for any combination of impurity 
concentration and potential. The AG curve plotted was 
obtained by fitting the unknown parameters to the crit- 
ical concentration for V = 20£, thus forcing the curve to 
go through the critical concentration calculated by the 
DCA. While for a given V and impurity concentration we 
cannot make a direct calculation of the pairbreaking pa- 
rameter entering the noninteracting AG theory and thus 
determine the critical concentration independently, the 



qualitative differences of our results from the AG curve 
shown arc obvious, particularly for small concentrations. 
The critical concentration calculated for strong impuri- 
ties agrees with the experimentally determined concen- 
tration for Cu-substitution by magnetic- and nonmag- 
netic impurities in LSCO[7[. 
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Figure 3: The dynamic spin susceptibility at Q = (0, 7r) for 
N c — 16, U — 4£ and V = 4f, at impurity concentration x = 
3%. The location of the peak, is a measure of the effective 
spin- wave exchange 2J„g [3a|. Inset: Spin coupling constant 
J as a function of V . 

Fig. [3] shows the magnetic structure factor S(Q,u>) at 
Q = (0, 71") of the system at temperature T = 0.087£ 
(rs 3T c ). In analogy with linear spin wave theory [38j, we 
note that the peak position of S(Q,uj) at Q = (0,7r) is 
a measure of the effective exchange coupling 2 J e g of the 
system. Therefore, we use S(Q = (0,7r),w) to extract 
J e jf of both ordered and disordered systems. We find 
that the rise of T c at low V is correlated with J e g of the 
system, as shown in the inset of Fig. [3] The initial rise of 
T c tracks the initial increase in Jeff- Then both T c and 
Jeff remain nearly constant up to V — t. Further increase 
in V causes J e g to remain roughly constant while T c 
shows a dramatic drop, indicating that the suppression 
at higher concentrations is indeed due to pairbreaking 
rather than pair weakening. 

For 3% impurity concentration, T c starts to drop at 
V ~ U/2, which is coincident with the formation of 
impurity-induced moments as shown in Fig. 2J A weak 
impurity does not induce any local moment in the sys- 
tem and increased J e ff causes a rise in T c ; increasing V 
causes formation of local moments and thereby enhanced 
pairbreaking. 

Discussion. Potential scattering due to weak local im- 
purities is expected to inhibit superconductivity, because 
the resulting isotropic scattering in momentum space 
causes <i-wave pair breaking and thus a reduction in the 
ii- wave order parameter. However, the suppression found 
here is in general weaker than expected from pairbreaking 
due to pointlike potential scatterers in a d-wavc system, 
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Figure 4: Squared local magnetization (see Eq. [3]) for N c 
— 16, and U = 4i, as a function of temperature for various 
impurity potentials and concentrations. The solid lines are 
guides to the eye. 



which follows the AG form. Our results suggest that 
the slowness of the initial T c suppression is due to the 
initial enhancement of the interaction by the impurities. 
These results are consistent with recent calculations J39J 
where the authors argue that an isolated impurity in a 
t—J model can enhance pairing locally. Since the instan- 
taneous part of t he p airing potential in the t-J model is 
proportional to </[40[, the local pairing and the transition 
temperature T c are enhanced. Here we have confirmed 
that for impurities, T c rises along with J. Intuitively, the 
increase in J can be understood by considering the 2 nd 
order exchange between two spins on sites with unequal 
energies, as discussed in Ref. [39J 

We are not aware of any experimental data indicating 
an actual increase or complete insensitivity of T c to in- 
creasing weak disorder in the cuprates, when doping is 
held fixed. It is not surprising, however, that our results 
overestimate the pairing enhancement effect of disorder, 
given the crude way in which disorder averaging has been 
performed here due to the current limitations on cluster 
size. Nevertheless, we regard these results as a strong in- 
dication that the observed slow initial suppression of T c 
in the cuprates, which has been remarked upon for many 
years, has its origin in large part in correlation effects. 
A point in the same general spirit was made within a 
different scheme for treating interactions in Ref. Il8i . 

Conclusions. We have studied the effect of pointlikc 
impurities in cuprates using the Dynamical Cluster ap- 
proximation. Our results show that for weak local im- 
purities, the superconducting critical temperature T c is 
weakly increased due to an average, impurity-induced 
enhancement of the antiferromagnctic exchange correla- 
tion J. With increasing impurity strength, local mo- 
ments start to form around the impurity site, causing 
more quasiparticle scattering, and the critical temper- 
ature plateaus and subsequently decreases due to pair- 



breaking in both potential and magnetic channels. The 
suppression of T c continues until the unitary scattering 
limit is reached, and T c remains constant. 

As a function of impurity concentration, T c is found to 
be enhanced by or insensitive to small amounts of disor- 
der, and although with large disorder T c is driven to zero, 
the suppression appears to be generally weaker than that 
predicted in Abrikosov-Gor'kov theory, where the slope 
of T c versus impurity concentration and potential is nega- 
tive for all concentrations and potentials larger than zero. 
Our results therefore strongly suggest that the observed 
slow suppression of T c is related to the strong correlations 
in the system neglected in the BCS approach to disorder 
in a d-wave superconductor. Together with the results 
of Garg et al.[25| and Andersen et al[26|], our work sug- 
gests a robustness of superconductivity in the presence of 
correlations against weak disorder in the charge channel. 
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